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ABSTRACT 

We set up a series of self-consistent TV -body simulations to investigate the fundamental plane 
of merger remnants of spiral galaxies. These last ones are obtained from a theoretical Tully- 
Fisher relation at z=\, assuming a constant mass-to-light ratio within the ACDM cosmogony. 
Using a Sersic growth curve and an orthogonal fitting method, we found that the fundamen- 
tal plane of our merger remnants is described by the relation R c 



1.48±0.01 r -0.75±0.01 
'0 J e 



which is in good agreement with that reported from the Sloan Digital Sky Survey R c 



1.49±0.05 



Re 



I e °- 75±0 01 . However, the i? 1 / 4 -profile leads to a fundamental plane given by 
In general, the correlation found in our merger remnants arises 



1.79±0.01 r -0.60±0.01 
T 



from homology breaking (V oc (7q, R g oc Kg) in combination with a mass scaling relation 
between the total and luminous mass, M oc . Considering an orthogonal fitting method, 
it is found that 1.74 < v < 1.79, 0.21 < 77 < 0.52 and 0.80 < 7 < 0.90 depending on the 
adopted profile (Sersic or i? 1 / 4 ). 
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1 INTRODUCTION 

The formation and evolution of elliptical galaxies continues to be 
an outstanding problem in astrophysics (e.g. Meza et al. 2003, 
Wright et al. 2003). In the current CDM models of hierarchi- 
cal structure formation, mergers are the common process through 
which larger objects form from small ones (Frenk et al. 1988, 
Kauffmann & White 1993, Lacey & Cole 1993). In this scenario, 
ellipticals are a 'natural' outcome (Toomre 1977) and accumulating 
observational evidence seems to support this (Schweizer 1998 and 
references therein). However, an important challenge to the merger 
hypothesis is to explain how this seemingly random merging pro- 
cess can lead to the tight correlations found in ellipticals (White 
1997). 

One of this correlations found in ellipticals is the so-called 
Fundamental Plane (FP); a linear relation, in the logarithmic space, 
among their effective radius, R e , their effective surface bright- 
ness, J e , and their on e-dimensional ce ntral velo city dispersion, 00 
iDiorgovski & Davislll987tlDressler et"al]ll987h : R c oc erjji*. The 
FP exhibits a general trend that follows from the application of the 
virial theorem to ellipticals: M oc V 2 R g , with M being the total 
mass, V 2 the 3-dimensional mean square velocity and R g the gravi- 
tational radius of the system (V 2 = 2K/M and R g = GM 2 /\W\, 
where K and W are the total kinetic and potential energy, respec- 
tively). Assuming that ellipticals are homologous systems in equi- 
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librium (V 2 oc a%, R g oc R c ) with a constant mass-to-light ratio, the 
virial theorem leads to 



R e oc a I e 



logi? c « 2 [log cr + 0.2/i c 



(1) 



where fi c = — 2.51og/ c . Deviations of an observed FP from the 
theoretical expectation are attributed to variations of the mass-to- 
light ratio with luminosity and/or to the breaking of homology in 
galaxies (e.g. Trujillo, Burkert & Bell 2004). 

The observed FP presents some differences depending on 
the wave- band spectrum and on the environment where galax- 
ies reside 1 l0reensen et al. ll996t|Pahre et alll998l : lMobasher et alJ 
ll99St iKelson et alJl2000l : bernardi et alj|2003l) . Recently, from a 
sample of about 9000 elliptical galaxies taken from the Sloan Digi- 
tal Sky Survey (SDSS), Bernardi et al. (2003) obtained a FP relation 
given by i? c (xai' 49±0 ' 05 /-°' 75±0 ' 01 ; where the J R 1 / 4 -profile was 
assumed. 

On the theoretical side, simulations of spiral mergers are nu- 
merous and have been extensively explored to study: the morphol- 
ogy of merger remnants (e.g. Mihos & Hemquist 1996), the line-of- 
sight velocity distribution (e.g. Bendo & Barnes 2000), and the ef- 
fect of the orientation of the angular momentum of discs (e.g. Naab 
& Burkert 2003). Other works have considered the effect of merg- 
ers on the properties of the FP (e.g. Capelato et al. 1995, Levine & 
Aguilar 1996, Bekki 1998, Dantas et al. 2003, Nipoti et al. 2003, 
Gonzalez-Garcfa & van Albada 2003). A common feature of pre- 
vious studies is that they have used models with properties chosen 
rather arbitrary resembling present day spiral galaxies. In partic- 
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ular, those works addressing the FP of merger remnants have not 
considered self-consistent models of disc galaxies. 

Disc galaxies also define a tight correlation between the max- 
imum of their rotation velocity V m and their luminosity, termed 
the Tully-Fisher (TF) relation, that can be written as L = AV£', 
where a and A are the 'slope' and the zero-point, respectively. The 
observed values for a lie between ~ 2.5-4, depending on the wave- 
band (e.g. Strauss & Willick 1995). In particular, the J-band TF 
relation is given by 



Mi - 5 log h = -21.00 - 7.68(log W - 2.5) 



(2) 



where Mi is the absolute magnitude of the disc, and W ~ 2 Vm 
is the 21 -cm hydrogen line width resulting in a slope of a « 3.1 
(Giovanelli et al. 1997). 

In this Letter we study whether merger remnants of self- 
consistent TV-body spirals satisfying the TF relation can lead to 
the FP of ellipticals. Spirals properties are obtained from the Press- 
Schechter formalism (Press & Schechter 1974) in combination with 
the model of disc galaxy formation of Mo, Mao & White (1998, 
hereafter MMW) and the TF relation under the ACDM model. The 
rest of this work has been structured as follows: in Sj2|we summarise 
the model used to set up the theoretical TF relation, the initial con- 
ditions and the orbital parameters for our simulations. In f|3] the 
results and a discussion are given. 



2 NUMERICAL MODEL 

In this section we summarise the method and properties of our mod- 
els to study the FP associated to merger remnants. 



2.1 Tully-Fisher relation of spirals 

To construct a theoretical TF we use the Press-Schechter frame- 
work of hierarchical clustering and the model of MMW for the for- 
mation of galaxy discs. In this model five parameters are required 
to obtain the radial scale-length of a disc formed inside a spherical 
dark halo. These are the circular velocity, V c , the spin parameter, A, 
the concentration of the halo, c, the fraction of disc to halo mass, 
nid and the fraction of angular momentum in the disc, jd- To ob- 
tain (V c , \,c,md, jd) we have proceeded as in Shen, Mo & Shu 
(2002). The ACDM cosmology was adopted to obtain the partic- 
ular properties of disc galaxies (with a mass density parameter of 
nn=0.3, a cosmological constant of Oa=0.7, the Hubble's constant 
h=0.7, and the perturbation power-spectrum normalization 0"8=1). 

We have used the theoretical TF at 2=1, an epoch thought to 
be close to the one where formation of the discs occurred (Pee- 
bles 1993), as our fiducial redshift to construct the numerical disc 
galaxies. This allowed us to consider remnants resulting from bi- 
nary encounters and to compare them with ellipticals that might 
have formed through these encounters. Present day massive ellipti- 
cals are probably the result of more than a binary merger, but their 
study is out of the scope of this work. 

We have selected only disc galaxies with a stability criterion 
e m ^ 0.9, where e m =Kn(GMd/7?d) _1//2 , and Vm is the maximum 
rotation velocity (Efstathiou, Lake & Negroponte 1982; Syer, Mao 
& Mo 1997). Galaxies were evolved in isolation during 2 Gyr, 
about twice the time needed for the first pericentric passage in the 
encounter. Despite the fact that our galaxy models match the previ- 
ous criterion some form a well defined bar (B) while in other cases 
they develop a transient bar (T) or remain stable (S) (see last col- 
umn of Table[j}. 
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Figure 1. Theoretical Tully-Fisher relation for an ensemble of spirals (dots) 
at 2=1. The dashed and solid lines correspond to the linear-square fit with a 
slope of a = —8.0 and to the observed relation (2), respectively. The prop- 
erties of discs of galaxies to be merged were selected from this ensemble. 
As an illustration, progenitors of mergers M01, M02 and A/03 of Table 1 
are shown. 



In Figure Q the TF relation, at redshift z=l, is shown. The 
straight solid line indicates the observed slope of the TF in the I- 
band of equation (2) at the present epoch. A sample of twenty ran- 
dom points from this TF relation were selected in order to set up 
the self-consistent TV-body disc galaxies. 



2.2 Galaxy models 

Our TV-body disc galaxy models consist of a disc and a spherical 
dark halo component. Following Hemquist ( 1 993), disc particle po- 
sitions are randomly drawn from the axisymmetric density profile: 

M d 



Pd(R,z) = 



■ exp(— R/Rd) sech {z/z d ) 



(3) 



being Rd and Zd the radial and vertical scale length of the disc, re- 
spectively; Zd was chosen randomly in the range (0.1 — 0.2)i?d- For 
the halo we have adopted a modified version of the NFW-profile 
(Navarro, Frenk & White 1997) with an exponential cut-off given 
by: 



Ph(r) = 



with 



4nr(r + r s ) ; 



exp 



(— + q 

\ T200 



(4) 



Oh 



exp( 



% /7rqexp(g 2 )Erfc(g) + | exp(q 2 )Ei(g 2 



being Erfc(x) the complementary error function and Ei (x) the ex- 
ponential integral, r s and V200 the scale and virial radii of the halo 
respectively, c = g _1 = (r s /r2oo) _1 the halo concentration, and 
Mh is the halo total mass. Finally, velocities are derived from Jeans 
equations. 

The parameters characterizing our disc galaxy models are 
listed in Table Q] where / d = M d /{M h + M d ). We have also in- 
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Table 1. Properties of merging galaxies. 



Merger Halo Properties Disc Properties 









A 


3 


c 


Ah 


fd 




2d 


Q 










[M Q ] 


[kpc] 












[kpc] 


[kpc] 




[kpc] 






M01 


4.02 x 10 11 


104.3 


0.063 


0.041 


7.63 


57126 


0.0498 


2.4 


0.39 


0.402 


5.8 


12000 


T 




1.25 x 10 12 


152.3 


0.056 


0.068 


3.84 


177571 


0.0586 


5.7 


0.47 


0.321 


13.8 


44267 


T 


M02 


1.39 x 10 11 


73.3 


0.021 


0.052 


8.94 


172403 


0.0312 


1.2 


0.12 


0.651 


3.0 


57436 


T 




6.45 x 10 10 


56.7 


0.028 


0.033 


6.61 


80000 


0.0236 


1.3 


0.13 


1.302 


3.1 


20000 


S 


M03 


8.65 x 10 10 


62.6 


0.078 


0.033 


11.12 


149138 


0.0402 


1.7 


0.20 


0.995 


4.1 


42093 


S 




4.64 x 10 10 


50.8 


0.031 


0.043 


12.38 


80000 


0.0357 


0.7 


0.12 


0.789 


1.8 


20000 


B 


M04 


2.36 x 10 11 


87.4 


0.036 


0.029 


9.17 


80000 


0.0393 


1.1 


0.19 


0.387 


2.6 


20000 


T 




1.75 x 10 11 


79.2 


0.048 


0.054 


5.29 


59428 


0.0372 


3.3 


0.58 


1.064 


8.0 


15671 


S 


M05 


5.42 x 10 10 


53.5 


0.053 


0.114 


15.95 


74829 


0.0817 


1.2 


0.19 


0.561 


2.8 


20000 


B 




6.02 x 10 10 


55.4 


0.034 


0.077 


12.05 


83124 


0.0541 


1.0 


0.14 


0.610 


2.3 


14245 


S 


M06 


7.49 x 10 10 


59.6 


0.098 


0.059 


7.91 


91546 


0.0638 


2.5 


0.35 


0.943 


6.1 


25000 


B 




6.68 x 10 10 


57.4 


0.078 


0.043 


13.32 


81742 


0.0420 


1.9 


0.23 


1.253 


4.5 


14327 


S 


M07 


4.88 x 10 10 


51.7 


0.074 


0.064 


11.02 


73643 


0.0635 


1.5 


0.25 


0.841 


3.6 


20000 


B 




9.48 x 10 10 


64.5 


0.047 


0.103 


11.79 


143062 


0.0656 


1.6 


0.31 


0.579 


3.9 


42119 


B 


M08 


9.77 x 10 10 


65.2 


0.032 


0.034 


6.56 


245462 


0.0239 


1.7 


0.20 


1.261 


4.1 


36000 


S 




1.02 x 10 11 


66.1 


0.023 


0.012 


7.80 


188124 


0.0127 


0.9 


0.12 


1.217 


2.2 


19725 


S 


M09 


8.11 x 10 10 


61.2 


0.099 


0.142 


10.87 


150000 


0.0818 


4.6 


0.73 


1.217 


11.1 


30000 


S 




8.33 x 10 10 


61.8 


0.122 


0.095 


10.01 


154050 


0.0791 


3.9 


0.56 


1.089 


9.6 


29686 


s 


M10 


4.74 x 10 11 


110.3 


0.112 


0.106 


11.39 


240000 


0.0807 


6.6 


1.29 


0.576 


15.9 


60000 


B 




6.99 x 10 10 


58.3 


0.071 


0.035 


9.94 


165992 


0.0396 


1.6 


0.19 


1.092 


3.9 


48689 


S 



eluded the radius Rq at which Toomre's Q parameter is normal- 
ized and the number of particles in the halo, A/h, and the disc, N^. 
Notice the range of concentrations of haloes and galaxy disc sizes. 

Our disc galaxy models do not include a central bulge com- 
ponent since, in the framework of MMW's galaxy formation, they 
tend to produce smaller disc sizes which require a better spatial res- 
olution and smaller timesteps and, hence, more demanding compu- 
tational resources. It is not clear how this central component will 
affect the FP of the remnants; it will be addressed in a future work. 



2.3 Initial conditions and experiments 

We have considered only parabolic encounters where the initial 
separation between galaxy centres is 7? s = 1.25(r2oo,i + ''200,2)- 
The pericentre for each encounter is randomly selected in the in- 
terval R p £ (5,20) kpc, consistent with the energies and eccentrici- 
ties of cosmological A/-body simulations (Navarro, Frenk & White 
1995). The relative orientation between galaxy spins is also taken 
randomly. 

To evolve the numerical simulations we use the parallel ver- 
sion of GADGET, an N-body/SPH code with individual timesteps 
(Springel, Yoshida & White 2001). Softening parameters of 
e c i=35pc and e h =350pc were used for the disc and the halo par- 
ticles, respectively. GADGET uses a spline kernel for the softening, 
so that the gravitational interaction between two particles is fully 
Newtonian for separations larger than twice the softening parame- 
ter (Power et al. 2003). 

The typical time of arrival at pericentre is about 1 Gyr, and 
all simulations were run for a total time of about 8 Gyr; this cor- 
responds to about the time span from z=\ to z=0 in a ACDM sce- 



nario. Remnants are already in a virial state by this time. All of 
our runs were performed on a cluster consisting of 32 Pentium pro- 
cessors running at 450 MHz (Velazquez & Aguilar 2003). Energy 
conservation was better than 0.25% in all cases. 



3 RESULTS AND DISCUSSION 

To analyse the merger remnants a constant mass-to-light ratio, T, is 
assumed; hence I c is proportional to the effective surface luminous- 
mass density E e . We first compute their 'photometric' proper- 
ties (R c , E c ) and then its central velocity dispersion <to inside an 
aperture of radius R c /8. The surface density profiles £(i?) were 
fitted using a Sersic profi le iS ersi 5 1 1 968t ICiotti & BertirJ IT999t 
Caon, Capaccioli & D'Onofrio 1993) and the R 1 ^ 4 profile, we also 
fitted the corresponding growth curves of the luminous part of the 
remnants: M L (R) = 2n f T,(R)RdR. 

To check how the photometric properties are affected by the 
adopted region to be fitted, we have considered an aperture with 
outer radius of 17.5 kpc as in Wright et al (2003) and for the inner 
radius of the fit, £, we use two different values: (1) the resolution of 
the luminous component in our simulations, £ ~ 2.8ed ~ 100 pc, 
and (2) £=300 pc. This also allows us to test the effects of the 
boundaries of the fitting region on the values of the exponents (a, b) 
of the FP (i? e oc (Tq^g)- Each remnant was 'observed' from 100 
random line-of-sights. A y 2 -minimisation by Levenberg-Marquart 
method I Pr ess et alll992t) was used to obtain R c and S e for each 
projection. The exponents (a, b) of the FP were fitted to the com- 
plete set of projected values of the remnants, using both orthogonal 
and direct fits. 

The resulting values for these exponents are plotted in Fig- 
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Table 2. Mean of projected values used to determine our FPf . 



Merger 






Parameters* 






(n) 




Wo) 


(Me) 






[kpc] 


[km/s] 




M01 


4.3 


9.0 


139 


-19.4 


M02 


2.2 


2.5 


75 


-19.4 


M03 


1.8 


2.2 


72 


-19.4 


M04 


2.6 


3.7 


95 


-19.4 


M05 


3.8 


1.7 


71 


-20.3 


M06 


2.3 


3.7 


67 


-18.8 


M07 


3.2 


2.5 


66 


-19.8 


M08 


2.8 


2.1 


45 


-19.1 


M09 


2.1 


6.1 


72 


-18.3 


AflO 


2.6 


5.2 


71 


-18.7 



* Here n is the index of the Sersic profile [ex exp(— ij 1 /")], 
and /z e = —2.5 log S c with E c in Mq /kpc 2 . 



ure|5|(Jeft panel). It can be seen how the values of parameters (a, 6) 
are sensitive to the fitting function (profile or growth curve), to the 
adopted profile (Sersic or i? 1//4 ), to the fitting method (direct or or- 
thogonal) and to the value of the radius f . The observed values in 
the r*— band of the SDSS have also been included (star symbols). 
In particular, the dotted lines refer to the exponents of the FP of 
Bernardi et al. (2003): R c oc cr 1 ,' 49 /" ' 75 . 

This last relation is quite well reproduced by the model with 
5=100 pc, a Sersic growth curve and an orthogonal fitting method 
(hereafter our fiducial model, denoted by FPf) which is expressed 
by the relation R c oc a i- 4 8±°-°i E -o.75±o.oi w ith a standard devia- 
tion orthogonal to the plane of 0.06. Uncertainties in the exponents 
were estimated using the bootstrap technique lEfron & Tibshiranl 
Il993l) . Table 2 lists the mean projected values of the quantities used 
to compute our FPf and the mean Sersic index. 

However, the exponents (a, b) obtained by fitting the merger 
remnants with the i? 1//4 -profile have problems matching the ob- 
served ones, and this is independent of the fitting method (or- 
thogonal or direct). For example, for £ = 100 pc and an orthog- 
onal fitting, the R 1//4 -profile leads to a FP described by R c oc 
a i.79±o.oi s -o.60±o.oi (which wi u be denoted by FP dv ), while for 

£ = 300 pc it is given by R c oc ct1' 65±0 - 01 E-°' 65±0 ' 01 ; notice that 
increasing £ tends to provide a slight better fit. An edge-view of 
our FPf for the merger remnants (filled squares) is shown in Fig- 
ure |2|(rigi3t panel). As a reference, the long-dashed line indicates 
the slope, av = 2, for virialised homologous systems with a con- 
stant mass-to-light ratio while the solid and dotted lines refer to the 
SDSS and to our FPf, respectively. Our merger remnants were ar- 
tificially displaced to coincide with the zero-point of the FP of the 
SDSS (Bernardi et al. 2003). 

The overall agreement of the parameter values (a, b) for the 
Sersic profile case with those of the SDSS relies on a combination 
of homology breaking (V 2 oc ctq, R g oc R]l ) and a mass scaling 
relation between the total and luminous mass as M oc M^. For 
instance, for our FPf non-homology is found to be V 2 oc o"o' 74 and 
R s oc R° 51 , with a mass scaling relation given by M oc Ml' 85 
and a negligible dependence of Sersic's index with the effective 
radius (n oc i?[?' 07 ), while for the FPav case, these relations are 
described by V 2 oc ct 1 ,' 77 , R e oc R^ 21 and M oc Ml' 80 . Note that 
a direct fitting method shows larger differences between our fitted 
FP and that of the SDSS (see Figure [3J-right panel). Despite that 
Sersic's profile fitting leads to a better match with observations; a 



consistent comparison will require an observational FP determined 
using this type of profile on the large data base of the SDSS. 

The above results suggest that, ignoring any contribution from 
breaking-homology relations, it does not appear feasible to scale 
the mass-to-light ratio, (M/L), as a simple power-law of the 
luminosity, L, (or equivalently of the mass M) alone; see also 
Bernardi et al. (2003). Perhaps, a more general relation of the form 
(M/L) oc RcO-qLI need to be explored. A more extensive study 
of the scaling relations of merger remnants, their non-homology, 
the distribution of luminous and dark matter, and their kinematics 
is under way and it will be presented in a future work. 

Our results support the idea that mergers of spirals, with prop- 
erties obtained from the Tully-Fisher relation within a ACDM cos- 
mogony, lead to the tight correlation expressed by the FP. 

In dissipational numerical simulations of mergers, it is found 
that a gas component can yield important differences in the stellar 
component of the remnants in comparison with their collisionless 
counterparts (e. g., Barnes & Hemquist 1996). For instance, a gas 
component deepens the potential well, affecting more the remnant 
kinematics than its luminous profile. However, it is not clear from 
these results how the FP is going to be altered; in particular, if its 
slope is going to be changed or just the dispersion of values around 
it. This remains an open issue. 
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Figure 2. (Left) Parameters (a, b) for the FP (R c oc uS/g) of our galaxy merger remnants. For comparison purposes, the values of the SDSS (star symbols) 
in the r*-band have been included. Horizontal dotted lines indicate the FP relation, JJ e oc Oq- 49 £~ ' 75 , quoted by Bernardi et al. (2003). To the left, these 
parameters were obtained using an orthogonal fitting method while to the right using a direct one. Here, S=Sersic profile, SQc=Sersic's growth curve, i? 1 / 4 =de 
Vaucouleurs profile, and Rqq its growth curve. Also, the effect introduced by two values of the inner radii § is illustrated. (Right) Edge-on view of the FP. 
Remnants are indicated by filled squares using the mean values listed in Table 2. The solid and dotted lines correspond to the SDSS FP (R c oc o-Q ,49 £,r ' 75 ) 
and to our FPf model. As a reference, the dashed line with slope ay = 2 corresponding to virialised homologous systems has been plotted. Remnant mergers 
have been artificially displaced to coincide with the 'zero-point' of the FP of the SDSS. 
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